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SU^IMARY 


TVis paper is s review of the nonequi libriura calculaf-an tecnniques 

for' finite rate catalyci.: acoai reco-abinacion on the taernal proce^tton 
naterials Th- ^e-hniq>'ies considered are the axisyoaetric viscous snock. layer 
methL -hrel-diraensional (3-D) reacting Euler equation solutions coupled uitn 
Ssiet^rrr-analog boundary layer method, and a recently Eevelopea conequili- 
briuu 3-D viscous shock layer method. 

The comparisons indicate a substantial influence of 
on the heating Co much ^ loier'" thln^ i/ the flow were in equilibrium 

°?thrLlefweJe’fullv catalytic. It is shown that all of the nethods agree with t 
measurements within about 10 to 30Z . depending on ^he location. 

and assumption about the catalytic -comoination [^/Tindward center- 

lations could predict the measurements uniformly ^ ^ calculations 

line for all fligat conditions. (Until now the 3-0 vis.ous .nock l.a>- - 

have only treated the noncatalytic wall.) 

X. • rhpt for a ziven flight condition the temperature measured on the 

It IS noted that tor a given -i-6 fjcrh flieht. The cause 

o?S in«^srirCrJL^!°bu^%rmay be dL to contamination of the surface, 
of this incr-ase -,talv'-i'-v or to a decrease in cmittance. 

causing an increase in catai/_i-y, or 

Nitrogen recombination was found to be significant early in the entry 
Nitrogen reco ^ s^ch as near -..ne nose. This 

atmosphere at higher velocities. 
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SYMBOLS 


ALT altitude 

C 3to!c mass fraction 

f heat flux adjustnent factors defined in eq. ^_2) 
hj total enthalpy 

k. Boltzmann constant 

catalytic recombination speed 
L length of vehicle 

T nass of atom 

P pressure 

q heat flux 

T temperature 

V velocity , ^ . 

VINE freestream velocity 

X axial distance from nose • 

Z geometric altitude 

Greek Symbols 

E emittance 

Y energy transfer catalytic combination coefficient 
Subscripts 

FC fully catalytic 

N nitrogen 

0 oxygen 

w wall 

® freestream 

ref reference condition or property 


INTimDUCTION 


The Space Shuttle orbiter is a hypersonic glide reentry vehicle that spends 
much of its entry time at relatively temuous altitudes in which chemical noneqioili* 
brium predominates in^the shock layer. Calculations have shown that both dissocia* 
tion ncnequilibrium^ and recombination nonequilibrium exist^ - The dissociated 
nonequilibrium exists in the inviscid layer and the recombination nonequilibri^um 
exists in the boundary layer. Verification of these phenomena has not been 
directly obtained; hcvcver, these phenoTsena are inferred by comparing heat transfer 
measurements with the reacting flowfield results. 

Although measurements of surface temperatures on the high temperature rettsable 
surface insulation (HRSI) tiles have be=:n made at numerous locations or. the 
orbiter, this paper only addresses measurements on or near the windwar*’. centerline 
of the lower fuselage because predictions of local flow conditions are much easier 
to obi-ain in this region. The presence of chemical nonequilibrium was made easier 
to verify because the HRSI tile glass coating (RCG) is relatively noncatalytic with 
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respect to ate® rocombir.ation and the associated dissociation energy acconmodatiot. 
Also of great importance in demonstrating the nonequilibrium flow behavior is the 
catalytic surface effects orbiter flight experiment of Stewart, Rakich, and 
Lanfranco^, whose initial results were reported in reference 4. Prior to the 
flight experiments , predictions of the noncatalytic nature were reported in 
references 1, 2, and 3 based on flowfield computations and arc jet experiments. 
Besides the results reported in reference 4 other calculations have been made ..or 
the RCG coated tiles and compared with flight measurements. Scott and Derry usee 
the reacting f lowfield/boundary layer method of reference 2 with measured energy 
transfer catalytic recombination coefficients of reference 6 and compared those 
predictions with flight measurements. Likewise, Shinn, Moss, and Simmoncs cc^utii 
heat fluxes 'ising an axisymmetric reacting viscous shock layer code with the recou- 
bination coefficients of reference 6 and showed better agreement with flight 
measurements. Recently Kin, Swaminathan and Lewis^ solved the 3-D viscous shock 
layer equations for the Shuttle geometry , ,ajd obtained encouraging results. 


This paper critically evaluates the various flowfield predictions, comparing 
the results of equilibrium and nonequilibrium flowfields coupled with reacting 
axisyrametric analog boundary layer solutions and the results of viscous shock layer 
solutions with flight temperature/heat flux measurements near the windward 
centerline for the Shuttle flights STS-2, 3, and 5. 


In the comparisons with flight heat flux measurements there is concern with 
two basic aspects of the predictions, the flowfield methodology and the surface 
catalytic recombination phenomena. The first aspect can be subdivided into 
dynamical and geometrical characteristics, and thermophysical properties and g^ 
phase chemical reaction kinetics. The second aspect can be subdivided into wall 
recombination rates of the basic thermal protection material, contamination issues, 
and knowledge gained from the catalytic surface effects experiment. All of tnese 
aspects are interrelated and the Shuttle flight does not provide an experiment in^ 
which each aspect can ’be controlled Independently. Numerical simulation is capable 
of single parameter variation, but confirmation of the results is difficult because 
of flight complexities and unknowns; particularly, there is no measurement o f t he 
chemical composition of the flow. This paper considers flowfield chemical ccwposi- 
tion effects (equilibrium vs nonequilibrium) , methods of solution (two- layer 
approaches and viscous shock layer approaches) , and surface catalytic recombina- 
tion rates, and it touches on possible contamination on the surface. The lsst>es of 
Incomplete chemical energy accommodation of catalytically-formed excited species 
and subsequent quenching are not explored. 


COMPUTATIONAL METHODS AND THEIR APPLICATION 


Five computational methods are considered here which are subdivided into 
applications of those methods, which are further subdivided into particular cases. 
These cas^s are sunmarized in Table 1 - 

The first two methods are axisynsaetric viscous shock layer methcxls of Moss^ 
and Miner and Lewis ^0. The next two are rwo-layer approaches . Rakich and 
Lanfranco^ treated the 3-D reacting inviscid flowfield and used the results 
reacting boundary layer edge conditions. Goodrich et al. solved the equilibrium 
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30 inviscid case and used their results as edge conditions for equilibriua 
boundary layer solutions. The fifth method is the 3-D nonequil ibrium viscous shoc/c 
layer method of Kim, Swamimathan and Lewis®. 

Shinn, Mess and Simnionds^ applied the Moss^ metaod with variable wall recotn— 
biration coefficients to Che Shuttle erbiter by approximating the Shuttle geometry 
with hyperboloids of revolution fitted by Zoby^^, They presented cases for various 
tires in the orbiter enfry and concluded the following. The Shuttle flight data 
indicates the shock layer flow is appreciably in nonequilibrium down to an altitude 
of 50 km. Scott's extrc.polated recombination rates® used in their viscous shock 
laT^r calculations result in good agreement with flight data forward, but not aft, 
on ihe vehicle. Better agreement aft is obtained if k^ = 100 cm/sec is used. The 
temperature of the surface during entry is 80 to 200 K less than if it were fully 
catalytic. 

Gupta, Moss, Siramonds and Shinn^® sioiiar-ly applied Che Moss^ method with 
various recombination coefficients and for a range of angle of attack of the 

orbiter. They found that a ±5^ variation in angle of attack does not affect the 

nonequi libriua heating appreciably at 75 and at 48 km altitudes. The temperature 

dep^dence of the oxygen recombination rate is not as steep as an extrapolation of 
Sccct's® data indicates. They concluded that a value of = 200 cm/sec seems to 

yield better agreement with the flight measurement of heat flux at certain 
locations and flight regimes. A 49T reduction in heating due to nonequilibrium 
effects was noted in the nose region at X/L = 0.025 and 7o km altitude. Nonequilib 
riun effects on the heating are not significant below about 65 km even though the 
flew may not be in equilibrium, xndicating that equilibrium boundary Li»yer methods ^ 

or heating correlations of the type suggested by Rakich et al- may be useful. 

The method of Rakich and Lanfranco^ was applied by Rakich, Stewart and 
Lacfranco^ to calibrate the results of an approximate method that uses equilibrium 
nomal shock isentropic boundary layer edge conditions in lieu of the reacting 
variable entropy edge conditions. This approximate method was then used to infer 
kyQ of the reaction cured glass (RCG)-coated high temperature reusable surface ' . « 

insulation (H3SI) tiles and to infer k^ of the iron-cobalc-chromia spinel (C742) 
coating used in the catalytic surface ^fects flight experiment tiles. They ; 

inferred that Islq = 80 cm/sec and assumed that = 0.3 ^wO for RCG at T„ of about 
llOQ K. Their ^talytic surface effects experiment . -^monstrated that the flow is ^ 

indeed in chenical nonequilibriua. Rakich' s method was also used by Scott^ with 
teirperature dependent recoabinacion coefficients inferred from arc jet measure- 
ments.®*^^ He used the reacting boundary layer code BLIMPK developed by Bartlett 
and Kendall and extended by Tong, Buckingham and Morse^®, This method resulted 
in higher heating than measured on the nose of the orbiter, but tended to predict 
or mderpredict the heating on the midbody. These results are presented here for 
comparison with other results - 

Reacting boundary layer calculations were made with equilibrium edge 
conditions provided by Gocjdrich et aJl.-*-^ along with different vail recombination 
assumptions. These results are presented here. 

Miner and LeWis-^ axisymmetric , reacting viscous shock-layer code was applied 
with various cataiytic wall assumptions and those results are likewise presented 
here . 
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The fifth ccuputat ional method considered herein for nonequilibrium flow 
calculations applied to the Space Shuttle was presented by Kim, Svaaiinathan and 
Lewis^. That recent paper addressed the windward side of the Space Shuttle, 
applying the 3-D nonequilibrium shock laye: method with noncatalytic boundary 
conditions. Their windward centerline results for two points in the STS-2 trajec- 
tory are presented here. 


MEASUREMENTS OF HEAT TuTZ TO SHUTTLE 


Surface temperature measurements of several instrumented HRSI tiles, distri 
bated along the lower surface of the orbiter, are considered in this paper. The 
flights considered are STS-2, 3, and to a limited extent STS-5 . Trajectory 
information was obtained from acceleration measurements on the orbiter and from 
atmosphere models calibrated by atmospheric soundings. The resulting best 
estimated trajectories (BET) were obtained frotn the Johnson Space Center, Hission 
Planning and Analysis Division. Heat fluxes were iaferrea from the measured 
temperatures by ccTmputing the corrected radiation equilibrium heat flux 

q = 1.06c f 

The factor 1.06 accounts for the fact that the cheniocouples lie about 0.38 tnm 
beneath the surface coating and for conduction in the tile. This factor was 
obtained from the aethod of Williams and Curry^^ who inferred heat tiuices from 
temperatures using an inverse thermal math model formulation*. Over the range of 
time in the trajectory and temperatures considered in this paper, a correction 
factor of 1 .06 is accurate to within 2 or 3 percent. 

When comparing the measurements of one flight ^th another or %rhen comparing 
calculations with measurements it is ^necessary to adjust the heat fluxes to account 
for differences in freestream conditions. Since the hypersonic stagnation poA.nt 
heating is approximately proportional to (p/pref^^^'^ ^ints were 

corrected by the ratio of that factor for the two freestream conditions, i.e., 

f =• ( Poo / Peref^^^^ (V»/Va,ref^^ 

The flight BET and the flight heating rates are used as reference cooditions when 
flight measurements are compared with calculations . The beat fluxes are then pre- 
sented in absolute units as obtained from equation (1). The factor probably 

accurate to within ± 3 % as verified by a comparison of calculations using the Miner 
and Lewis^® code. All the comparisons were made for an angle of attack of about 40 
degrees . 

To determine the consistency from flight to flight the bottom centerline heat 
flux measurements for STS-2, 3, and 5 are compared at three different times in the 
flights as shown in Figures 1, 2, and 3, respectively. The corresponding free- 
stream conditions are given in Table 2. It is seen that the flight— to— flight 
repeatability Is about 15-30Z and the standard deviation about the mean at each Z/L 


* The author is grateful to S. D. Williams of Lockheed Engineering and Management 
Services, Co. for calculating the heat flux for this determination. 



r^°?‘ neat fluxes are consistently lover than the other two 

he almost all locations. Tne reason for this is not understood, but i: ,ay 

3 ,-v, ^ ° ^ c ange iti catalycity resulting froa contamination of the surface or 

‘han emittance. The heating rat? correction for these cases is r.o larger 

^able 2. The measurement at X/L = 0,695 set-mt- anomalously low 
ere ore, it may be a bad measuremeat . The rather large discrepancy ia t-ie 

”d~;nod " “‘S'-" « 


SHUTTLE CENTERLINE PRESSURES 


Pressure measurements during the time of high heating on the orbiter were 
nea onlv rfurino ctc— i c -m . . « 


, . 1 ^ . — UL uigii neating on cne o 

ne only during STS-3 and 5. These measurements normalized by 

Dresen^pfi in ■FTnn**cv f, ^ .. 


^ , . ' , ^ mcaauicmenLS nonaaiizea DV V— aro 

Jhirtho calculatid -by three methodl. It is seen 

® repeatability is very good at almost all locations. The 
almo!r ^ decreases very rapidly m the forward lOZ of the vehicle then remains 

flowfield rr I/L of 0.1 to 0.4, rising slightly at X/L - 0.8. The 3-D 

flowfield calculations of Rakich and LanfrancoZ and Goodrich et al.^1 agree witi the 

the viner ard measurements. Likewise, the calculations using 

the ...ner and Lewis ^0 code agree within about 5Z except at X/L =■ 0.1 where the 

ar::::rL:ntrr'" discrepancy at X/L ^ 0.1 ma^ «s“t1rom 

ario^an" s" I bias of unknown amountlS. The instrument- 

readinvs^ ^ processing of the pressure measurements only result in positive 

. " exi8cca._s of a negative bias was indicated by a measurement that did 

caicula-inn« ® associated with the negative bias is saall, then the 

Jhe ° Although 3 diract comparison of 

tte TloTflTlTc T T " possible that the geometry descriptions ia 

flowfield codes do not adequately describe the vehicle as actually built- 

°ve^"'dica,'r,L t»e „pid „pa„,io„ „Lnd (^4 

-verpredicting the pressure (ana heat flux) near X/L =0.1. 


HEAT FLUX COMPARISONS 


Ir. the following comparisons the author has used the results of others and in 
.p "■ Of Pth-rs to „k, ptopopt cdlcolotiopp. 

author IS responsible for auy error or misapplication of the method, not 
the developers of the nethods. 


Equilibrium and Nonequilibrium - Two Layer Methods 


calr. shown in the past that the heat flux predicted by equilibrium 

maJeW ^ calculations with a fully catalytic wall are approxi- 

y qual. To verity this for the two boundary layer methods cosidered here 
comparisons made between equilibrium results of G^drich et al.H and'thr 

IrrlJ* n Lanfranco2 method. The freestream conditions 

with two nnnp ^ f^^gure 5 the Goodrich equilibrium prediction is ccmparec 

with two nonequilibriua boundary layer cases with fully catalytic walls. (Fully 
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catalytic here means Yq ^ = 1, vis-i-vis ky * or Cq^ * = 0.) I-^*e edge 

conditions in the latter two cases are Goodrich equilibriun and Rakich and 
Lanfranco nonequi librium . Given the same edge conditions it is seen "hat the 
equilibrium boundary layer calculation is about 15Z lower than the all-nonequi li- 
brium calculation over the entire length of the vehicle. Evidently, the transport 
of chemical en’^rgy by diffusion is more efficient in this case than via conversion 
of chemical energy to therral energy hicn is then transported to the wall via 
conduction. The opposite result was obtair.e'’ by Shinn et al.' who found the equilib- 
rium viscous shock layer resulted in higher heating than the equilibrium catalytic 
wall nonequil ibriun case. 


The comparison of the reacting boundary layer with equilibrium edge conditions 
’^ersus reacting edge conditions indicates that on the nose the’^e is very little 
difference between the two cases, whereas on the midfuselage the nonequilibrium 
edge conditions results in about 151 lower heating. The nonequilibrium edge case 
with a fully catalytic wall is very close to the all equli^irium calculattcn aft of 
X/L = 0.2. The latter agreement does not stem from the flow approaching 
equilibriun downstream because the equilibrium nitrogen atom concentration both at 
the edge and in boundary layer is greater than the nonequilibr iuta concentration by 
a factor of abcjut 1.3 in this case. Moreover, it was shown in reference 1 that the 
boundary layer is virtually frozen. 

A conrarison of the sxi symmetric reacting viscous shock layer method of Miner 
and Levis and the reacting two -layer approach of Rakich and Lanfranco-^ is made in 
figure 6 where the nonequilibrium boundary layer result for a fully catalytic wall 
lies above the viscous shock layer results by about 20Z on the nose. Agreement 
improves to within about lit at X/L = 0.55. The effect of edge conditions is about 
lot cr less everywhere for a fully catalytic wall. The results for a ncncatalycic 
wall are given in figure 7 where it is seen that the boundary layer heat flu:: :s 
about 30% greater than the viscous shock layer heating on the nose, but imprcves to 
about 10% at X/L = 0.55. Agreement of the reacting viscous shock layer resuxts and 
the reacting boundary layer with equilibrium edge conditions is within ab’ut 10% 
everywhere along the body. The equilibrium edge condition results fall below the 
noriequil ibriucn results on the nose, but they are very close farther aft. 


Attention is now turned to a comparison of the axisymnetric viscous shock 
layer method of Moss^ as applied by Shinn et al.^, and the two-layer method of 
Rakich and Lanfranco^ applied here for a lower velocity and altitude situation. 

It is seen in figure 8a that the axisymmetric nonequilibrium viscous shock-layer 
with equilibrium catalytic vsll (ECW) and the equilibrium ascous shock layer agree 
quite well (within about 5%). They also agree quite veil with the Goodrich^^ 
equilibrium two-layer result. It is seen that the fully catalytic ronequilibrium 
two-layer results are greater by about 10-20% which is the same as noted for case 
1 . Agreement in figure 8b for the noncatalytic case is worse than the two-layer 
results, being about 20-40% higher than the axisynmetric viscous shock layer 
results of Shinn et al.^ and present results using the Miner and Lewis^^ cooe. The 
latter results seem to indicate that heat transfer by atom diffusion is mere 
important in the viscous shock layer. This is consistent with the somewhat higner 
degree of dissociation, especially the nitrogen, associated with the viscous shock 
layer calculation. The reason for the differen*ces in atom fraction in the two 
methods is not understood since the reaction rates used in both methods were 
essentiaHy the same. 
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The recent 3-D nenequi 1 ibrium viscous shock layer results with a noncatalycic 
wan arc also giv--n in figure 8b. The heat flux at X/I = 0.1 is closer to the 
axisyrmetric VSL results, but <ioes not decrease as rapidly townstreaa. In fajt, 
tn^ 3-D results are hightr than tSie tvo-layer results afr of ll/L = 0.2. The 
chenical reacti'.yn tiodel-^ in all the ^aeLhods is virtually the sacie, (except ions 
are neglected in the Rakich and Lanfranco oiethcd). Therefore, the differences seen 
are dost likely due to differences in the computat ic**;al caethod or the geometry. 

Comparison of Measured and Caicuy^^e^ Hect Flux 

Attention is now turned toward a coapariton of the calculated htot flux and 
the measured values along the lower surface centerline The comparison is two 
times in each of two flights. The particular times were s elected r. > match the 
velocity and density as closely as practicaole to the conditions used in the bound- 
ary layer predictions for cases 460 and 65.1 in Table 3. It vas not possible to 
simultaneously match both velocity and density. The resul t\ng 'heat fluxjs were 
adjusted for the mismatch by the factor f of equation (2). As mentioned sarlie*^ 
the measured heat fluxes were inferred from the measured temperatures using 
equation (1) where E ~ 0.85. 

Several choices of catalytic recombination coefficients were used as wall 
boundary conditions for the two layer and the axisycrmetric shock layer calcula- 
tions. The energy transfer catalytic recombi nation coefficients for nitrogen and 
oxygen recombination on the RCG tile coating are presented I- figure^ 9 and 10, 
respectively. The coefficients presented are those found ia references 6, 7, and 
12. In those cases where a catalytic speed k^, was given the recombination 
v,oe *. f i c i ent is plotted as a dashed line, the lengtn of which indicates the tempera- 
ture range over which k,. was used, wiiere 


^ 2 r ci/kT 


It is seen in figure 9 that the inferred values of of reference 4 are 1 to 
0.2 times the values of reference 6. This lack of agreement is not surprising 
since k*/^ was a s sumed to be 0.3 times In reference 4, The values of (see 
figure iO) of references 4 and 6 a^.ree within experimental accuracy at the higher 
temperature range. At lo r temperatures = 80 cm/sec is about a factor of 
1.5 to 3 higher than the extrapolation of reference 6, depending on temperature. 

Extrapolating Co such a low temperature could be inaccurate, but the extra- 
polation is generally consistent with other recombination measurements (see figure 
6 of reference 4). Since the temperatures measured or* the Shuttle fell mostl-" in 
the lower temperature range 900-1100 fC, the predictions of heat flux using the 
k^ = 80 cm/ sec would result in higher heating except on the nose or earlier in 
time v;hei.e the nrtrogen carries a larger part of the dissociation energy. 

A comparison is made in figures 11-14 between the measurements and several 
calculations for STS-2 and 3 at iwo times in the trajectories. The measurements 
are near the bottom centerline of the vehicle except for a fev points that are ibout 
1.3 m off the centerline. Tb.is comparison between the calcuja * ' ->ns lud the measure- 
ments is typical for all times and both flights. In the higrie. :;titude cases 
(figures 11 and 13) the viscous shock layer methods with the temneranire depenc‘=‘n^ 

6 cf , . - - . 


values 
with k 


and 




yield better agreement for X/L < 0.3. The 


wO 


= 80 cm/sec also agrees with the measurements at X/L 


;o layer method 
. 5 . \t the lower 
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altitude (figur6S 12 and 14) the two“layer methods yield better agreement for 
X/L >0.2. 

In the higher altitude cases (figures 3.1 and 13). *:ae “onequi librium 
axisymnetric viscous shock layer methods with the te^;*rrac-..re dependent Yq and 
yield good agreement at X/L<0.3 and the nonequilibrium two layer method with 
k^ = 80 cm/sec and k^^ * 24 cm/sec yields good agreemenr for X/L <0.5. For the 
lower altitude case this two-lay^^r approach yields better agreement for X|L>0.2 
than the axisymmetric viscous shock layer methods. The nonequilibrium two-layer 
method using temperature dependent Yq and results are about 30% higher than the 
measurements on the nose area for all cases presented here, but the agreement 
improves toward the mid-vehicle and at i.,vcr altitude. It is apparent ’“hat the 
two-layer approach predicts higher heat fluxes for given wall boundary conditions 
than the viscous shock layer approaches. This may be due ^p^rt to the VSL having 
a slightly higher level of dissociation as well as to differences in the flowfield 
dynamics . 

In comparing the 3-0 nonequilibrium calculations of Kim, Swaminathan and 
Lewis^ with other noncatalytic predictions, one sees that the heat flux does not 
decrease as rapidly down the vehicle as do the axisymmetric viscous shock layer 
calculations and the two layer calculations. This indicates a possible influence 
of geometry and cross flow that is more adequately accounted for in the 3-D viscous 
shock layer model. In figures 11 and 12 the 3-D viscous shock layer results of 
reference 8 tend toward better agreement with the measurements than the other 
calculations aft of X/L = 0.6. This 3-D approach should be further investigated 
with appropriate cinite rate recombinaticn coe f ficients . 

The incx'ease in measured heat flux above the calculations on aft half of the 
vehicle and especially for the later flight may also have ocher explanations. The 
increase could be due to increasing recombination rates, but that would be incon- 
sistent with the mearurements on the forward part of the vehicle unless the^^ft is 
contaminated with a catalytic material. This is possible because the adhesive used 
Co bond tiles to the structure contains various metal oxides, particularly iron 
oxide which is known to be r.ighly catalytic. 

The two— layer methods have been used to calibrate faster and more flexible 
codes to provide heat fluxes and other properties over a wider range of conditions 
than for which the two-layer methods were applied. The nonequi librium results of 
Rakich and Lanfranco^ have been used by Rakich, Stewart and Lanfranco and by ScouL 
and Derry^ . One of the weaknesses of these applications is the inability to 
properly account for variations in the flowfield chemical composition as parameters 
fuch as the angle of attack, freestream dcifsity and velocity differ from the few. 
cases available from the 3-D Euler solutions. The axisymmetric shock layer codes 
have tne advantage that they are mere flexible in running cases because of their fast 
computaticn rime. Oupta et al.^^ investigated the influence of smalj. variations in 
angle of attcck on the nonequilibrium heating and iound the influence on heat flux 
to be small . A larger percentage variation was observed for lower velocities. 

This may be d\ie to greater temperature’ sensitivity of the level of oxygen dissocia-. 
tion at lower temperatures rising in the lower altitude case. Small changes i.n the 
component cf velocity normal to the shock wave associated with the change in angle 
of attack result in temperature changes in a range in which the oxygen dissociation 
is very sensitive. However, since the general sensitivity of absolute heat flux to 
angle of attack is small, the approximations made in references 4 and 17 should not 
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be very significant in this regard. 


The disadvantage of the axi symmetric viscous shock layer methods is that they 
are only capable of handling bodies of revolution and they rely on angle of attack 
simulation via changing the body profile. Cross flow or transverse body curvature 
IS therefore quite limited. Fortunately, for the present work this has not been a 
strong limitation but it may explain why there is disagreement with measurements on 
the art of the vehicle. 

The 3-D viscous shock layer approach does away with those approximations, but 
suffers from the requirement of a shock shape as input (as do the axisymmetric 
viscous shock layer solutions). The 3-D version has only recently been developed 

and will require further work to compare with measurements before its adequacy will 
be known . . » . 


The 3-D inviscid solution method coupled with boundary layer solutions 
requires very much computer time to obtain the inviscid flowfield and requires 
assumptions about how far into the inviscid flow from the body to go to obtain 
boundary layer edge conditions. Choosing the boundary layer edge Loo far into the 
lower entropy flowfield will result in heating predictions that are too high. This 
may be the reason that the two-layer methods predicted higher results than the 
axisymmetric viscous shook layer method for the noncatalytic case at X/L>0,02 and 
for the fully catalytic case at 0 ,02 < X/L < 0 . 2 . The noncatalytic case is more 
sensitive to the dissociation level which is higher in the flow from the normal 
shock region. 

catalytic reconbinat ion rates from the flight measurements is made 
1 iculc .or several reasons. First, the flowfield is composed of oxygen and 
nitrogen atorab in varying amounts according to tbs vehicle trajectory end location 
on the vehicle. If one chooses a lower velocity condition where very ’ittle 
nitrogen is dissociated then it may be possible to infer k^. However, we have 
seen a Oight-to-f light measurement uncertainty of at least 15% and prediction-to- 
prediction variation of the same magnitude. Heating uncertainties of this 
Mgnitude result in uncertainties on the order of a factor of 5 (s'ee reference 
131. Therefore, such a procedure should be used with great caution. This 
illustrates the need for careful ground experiments or great fidelity in the flight 
heat flux calculations to obtain accurate recombination coefficients. The ground 
measurements of Scott • , as with any ground measurements, require either precise 

heat flux calculations and/or a reliable reference surface with which to compare 
the heating. Even then accurate results are difficult. Fortunately only 
moderately accurate recombination coefficients are required to calculate reasonably 
accurate heating rates. 

To ascertain the nitrogen recombination coefficients from flight measurements 
is almost impossible without knowledge of the coefficients for oxygen because the 
oxygen atom is always an important species in the flow whenever there is any nitro- 
gen dissociated. If the flowfield and k^ were known accuiately as a function of 
T,^ then it might be possible to infer k^j^. 

Heating to Highly Catalytic Tiles 

Attention is now turned to the results of the Ames Research Center*s catalytic 
surface effects orbiter experiment^ >4. Not only did this very significant 
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experiment demonstrate the noncatalytic nature of the RCG coated tiles, but it also 
mav give some clues as to the variation of dissociation in the boundary layer an 
the wall recombination rates. On STS-2 two tiles were painted with a highly cata- 
lytic material, iron-cobalt-chromia spinel (C742) , developed by Stewart et al. ^at 
the Ames Research Center. The predicted and measured heat fluxes in the vicinity 
of the rvo C752-coated tiles on the bottom centerline of the orbiter during STS-- 
reentry are given in figures 15 and 16. The measurements were obtained at 
475 sec after 122 km altitude was reached. At the forward location near the nose, ^ 
the two— 1 aver calculation using the method and recombination rates of Rakich e- a 
yields Che best agreement with the measurements. This should be the case because 
the recombination coefficients were inferred from the measurements at this location 
and approximate entry time. Also shown is the same calculation but using 
recombination coefficients obtained from arc ^6t measurements » . As see in 

figure 15, (tne forward location), the increase in heat flux on the C742**coated 
tile is larger for the Rakich recombination rates than for the recombination rates 
of referer.cer 6 and 14, even though the latter rates for C742 are larger. The 
reason for :his behavior is that, due to the higher RCG recombination rates of 
reference 6, the boundary layer is depleted of atomic nitrogen and oxygen so that 
when the flow reaches the C742-coated tile there is not as much chemical energy 
available for transfer to the highly catalytic tile. A similar behavior is seen at 
X/L == 0.4 in figure 15. Since the recombination rates of references 6 and 14^ 
increase with temperature, the upstream edge of the C742-coated tile sees a high 
beat flux that decreases rapidly because of depletion of atoms in the boundary 
layer and thia lead to further reduction in recombination rate along the tile as 
the temperature decreases. 

The agreement between the axisymmet ric viscous shock layer method and the 
boundary layer method is not very good on the C742-coated tiles. The heat 
drops much more rapidly, possibly because of more rapid depletion of atoms in the 
boundary layer, than the boundary layer method predicts. There also seems, to oe 
some sensitivity of the heat flux distribution along the tile to the streamwise 
nodal spacing used in the calculation. 


CONCLUDING REMARKS 


This paper has attempted to evaluate the current state-of-tne-art nonequili- 
brium flow tools applied to the Space Shuttle. From this discussion the importance 
of nouequilibriuni phenomena to the Space Shuttle reentry heating has been assessed. 
Since the inception of the design of the Space Shuttle over fourteen years ago 
there have been developments in the heat flux prediction methodologies. Initially 
nonequilibrium and surface catalysis effects were ignored. Ihis led to a desip^^ 
that exceeded the requirements in many areas, but also resulted in an added margin 
pf safety in other areas that proved beneficial. 

It was found that the heat fluxes Taeasured on Che windward centerline of the 
orbiter tended to increase from flight to flight. Roughly, a 20^ cnange WaS noted 
from STS-2 to 5 at most of the thermocouple locations, indicating changes in 
surface properties such as emittance or catalycity. 

The nonequilibrium heat flux methods that have been developed and the 
catalycity measurements obtained over the past decade have improved the prediction 



capability from a 20 to 1002 overpredictioii for an assumed fully catalytic surface 
material to an accuracy of about 10 to 302 for nearly aoncatalytic materials, e.g. , 
the RCG coating on HRSI . These methods are the two-layer inviscid 5-D reacting 
flowfield coupled with the reacting boundary layer, and the reacting viscous shock 
layer solutions. The application of these methods may result in less reliance on 
wind tunnel measurements which cannot simulate the high enthalpy reacting flows 
associated with orbital reentry. Indeed the calculations are necessary for such a 
s imu lat ion . 

As the comparisons of the predictions with the measurements from the Space 
Shuttle flight tests have shown, we are now in a position of refining the predic- 
tion techniques and determining those phenomena that will be of significance for 
the design of future reentry vehicles such as an aerobraking orbital ^t^an^fer 
vehicle (AOTV) . 

Although nonequilibrium calculation te>-haiques using finite rate catalycity 
wall boundary conditions has significantly improved the prediction capability, none 
of the methods yields good agreement uniformly for all locations and freestream 
conditions. This points to the need for further refinement in these methods. The 
3-D viscoua approaches in particular should be pursued since the trends of the 
heating profiles tend to be better than for the other methods. 

Nit.'cgcn recombination is seen to be a very important phenomenon, particularly 
on the nose and elsewhere at the higher velocities. This means that the accuracy 
of the nitreg^’n recombination coefficients is important to the heat flux predic- 
tions in those 3reas. Since the AOTV enters the atmosphere at higher speeds and 
remains at highe'* altitudes where nonequilibr iuia flow dominates, the nitrogen gas 
and surface reactions will be especially important. 
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TABLE 

1 - METHODS AOT APPLICATIONS PRESENTED 

IN PLOTS 


Application 

No. 

Method 

Application 

Wall 

Boundary 

Condition 

Flow Field 

Boundary 

Edge 

Condition 

Layer 

Cheo' stry 
of 3.L. 

1 

Moss^ Shinn, eC al^ 

Scott^ 

VSL+ Hooeq. 

nTa 

N/A 

2 

11 

11 

ECW* 

II 

II 

11 

3 

11 

It 

Noncata . 

If 

II 

It 

4 

11 

ir 

Equilib . 

VSL Equilib. 

ir 

II 

5 

Miner & 
Lewis^O 

Present 

Scott^ 

VSL Reacting 

II 

It 

6 

II 

11 

Fully Cata 

II 

II 

It 

7 

II 

II 

Noncat^. 

II 

If 

It 

8 

Goodrich , 
eC al^^ 

Goodrich , 
et ai^l 

Equilib . 

3-D Inviscid 
Equilib . 

Equilib . 

Equilib . 

9 

Rakich . 
et al^ 

Present 

Fully Cata 

. 3-D Inviscid 
Noneq. 

Noneq , 

Noneq . 

10 

IT 

II 

Scott^ 

II 

It 

II 

11 

• 

II 

n 

Koncata. 

II 


II 

12 

Goodrich , 
et al^^ 

n 

Fully cata 

II 

Equilib. 

M 

13 

II 

*1 

Scott^ 

II * 

If 

II 

14 

11 

m 

Noncata. 

II 

II 

II 

15 

Pakich j 
eC al^ 

11 

Rakich . 
et al^ 

If 

Noneq . 

II 

16 

Kim, et al^^ Kim,et al 

Noncata. 

3-D VSL 

h/a 

N/A 


Noneq . 


* ECW = Equilibrium Catalytic Wall 
+ VSL * Viscous Shock-Layer 



CJ?iCiTv',*L F3 

OF POOR QUALITY 

TABLE 2 - FREESTREAM CONDITIONS*^ FOR CORRESPONDING 
HEAT FLUX MEASUREMENTS 



t ioc^ 

V, 


a 


z 



Flight 

sec 

kc/sec 


deg 

' 

ka 

.Nj’kg 

kPa 

STS- 2 

i75 

7.16 

.412.4 

40-37 

1,000 

74.7 

25 .6 

, , , 


700 

6-57 

.807.4 

:).99 

1.000 

70.2 

21.5 

3.49 


lO'O 

4.55 

.402-3 

40.56 

1, 000 

57 .2 

12. 4 

3.35 

STS- 3 

400 

7.29 

.394-4 

40.02 

1 .033 

75.1 

26 .6 

2.09 


700 

5.29 

.113-3 

39.58 

1 .037 

68,2 

1 ? . 8 

4.46 


960 

4.5S 

.417-3 

40.72 

1 .032 

57.6 

1 J.5 

8.75 

STS- 5 

400 

7.17 

.408-4 

40.05 

.998 

74.9 

25-7 

2.09 


700 

6.19 

.lo4-3 

40-71 

.950 

61.9 

19.2 

3.99 


950 

4.56 

.441-3 

39.19 

1 .O-'S 

54.0 

ID. 4 

9.16 

a f - 


(v/v „ 

f)3 

Factor used to 

adjust heat flux 

relative 


STS~2 condition « based on stagoacion poiac 
theory. 


b Tiae from entry interface CZ ■ 122 ka) 

c Best Estimated Trajectory 


TABLE 3 - FREESTREAM CONDITIONS FOR CALCULATIONS 


Case No. 

1 

1 

2 

460^ 

650* 

Method 

8L 

VSL 

BL 

VSu 

VSL 

Velocity, ka/ 5 

7.62 

7.62 

6.514 

7.20 

6,73 

Attitude, ka 

75.0 

75.0 

68-9 

75. o; 

71.3 

Angle-of-Attack , deg 

41 .4 

41 .4 

40.2 

40.0 

39.4 

Density, 

3.795-5 

3.974- 

5 9.28-5 

3-81-5 

6.33- 

Total Fnthalpy, Mj/kg 

29.0 

29-0 

21 .3 

25.9 

22.6 

Temperature, It 

197. 

197. 

221 . 

193. 

205. 

Stagnation Point Pressure. kPa 

2.20 

2.3* 

4.10 

1.93 

3.09 

N' ^e Radius , m 

0.814 

1 .342 

0.814 

1 .276 

1-233 

Hyperboloid angle^, deg 

- 

4-. 2 

- 

40-75 

40-20 

Lewis No. in Shock Layer 
Present Calculations 


1 .0 


1 -0 

1 .0 

Reference 11 

- 

1-4 


1 .4 

1 .4 

Reference 16 

- 

1-4 

- 

1 .4 



a These freestrcan condirions are tr.e sane as the one in reference ' 

for STS-2 times correspcF.diog to the case niniber and the same as ..aaes 2 and 3 
in reference 8. 

b Not applicable to 3-D VSL. > 
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Figure 7.- Comparison of nonequilibrium axisyimnetric 
viscous shock layer with two-layer calculations 
having nonequilibrium and equilibrium edge con- 
ditions for noncatalytic wall. 
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Figure 10.- Fuergy transfer catalytic recombination coefficient 
for oxygen on RC(l-coate<i tiitSI. 
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Figure 11.- Comparison of calculated and measured Shuttle heat 
fluxes for STS-2. t = 475 sec; altitude - 74.7 km. 
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of calculated and measured Shuttle, heat 
t = 475 sec; altitude - 71,29 km. 
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Figure 13.- 


Comparison of calculated and measured Sh ittle heat 
fluxes for STS-3. t * 400 sec. 
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Figure 14.- Comparison of calculated and measured Shuttle heat 
fluxes for STS-3. t = 700 sec. 
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Figure 15.- Heat fluxes on C742-coated tile at X/L = 0.163 
for STS-2. t = 650 sec. 
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Figure 16.- Heat fluxes on C742-coated tile at X/L = 0.40 
for STS-2, t = 650 sec. 
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